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Abstract Using an effective one-band Hubbard model with disorder, we con- 
sider magnetic states of the correlated oxide interfaces, where effective hole 
self-doping and a magnetially ordered state emerge due to electronic and ionic 
reconstructions. By employing the coherent potential approximation, we ana- 
lyze the effect of random oxygen vacancies on the two-dimensional magnetism. 
We find that the random vacancies enhance the ferromagnetically ordered state 
and stabilize a robust magnetization above a critical vacancy concentration of 
about c = 0.1. In the strong-correlated regime, we also obtain a nonmonotonic 
increase of the magnetization upon an increase of vacancy concentration and a 
substantial increase of the magnetic moments, which can be realized at oxygen 
reduced high-T c cuprate interfaces. 

Keywords Oxide interfaces • Titanates • Magnetism ■ Interface reconstruc- 
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1 Introduction 

In complex thin-film transition-metal oxide heterostructures, interface phe- 
nomena can substantially affect the electronic properties and lead new elec- 
tronic ordered states not observed in the bulk constituents 
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Prominent examples are the titanate interfaces in the heterostructures LaAlOa/SrTiOs, 
LaTi03/SrTiC>3 and similar compounds [8,9 . These structures have recently 
attracted much attention which is focused on the discovered two-dimensional 
electron liquid state stabilized by the interface electronic reconstruction. The 
magnetoresistance and torque magnetometry measurements [TUIITTIIT^] demon- 
strate a hysteretic behavior in the field-dependences of magnetoresistance and 
superparamagnetic behavior which suggests the existence of ferromagnetic 
puddles in the samples of LaA10 3 (LAO) grown on SrTi0 3 (STO). 

The observation of the magnetism at the LAO/STO interface triggered 
an intense exploration of the role of impurities in the formation of the mag- 
netically ordered state at these interfaces [S l fTU t fT ^ rTT irT^irT^lfTKlfT^i rTT ] . In par- 
ticular, recent first-principle studies demonstrate a vacancy-related magnetic 
exchange splitting of 3d states of interface Ti atoms [TBlfTT] . In this case, the 
conducting electrons, which emerge due to the charge compensation of the 
polar discontinuity occupy d xy bands, whereas the vacancy-released electrons 
occupy e g states shifted below the Fermi level due to the interface orbital 
reconstruction [161117] . STM, cathode luminescence studies and conductivity 
measurements provide strong support for the existence of the oxygen vacancies 
in STO layers of LAO/STO heterostructures [THUS]. 

As the DFT calculations are unable to access a realistic random low- 
concentration distribution of impurities, we consider here an effective two- 
dimensional one-band Hubbard model on a square lattice with disordered ran- 
dom vacancies. In this model, each vacancy introduces an exchange splitting 
of the local 3d xy state of neighbouring Ti atoms, in this way stabilizing a 
ferromagnetic order through the electronic transfer term. 



2 One-orbital model of random oxygen vacancies 

The two-dimensional electronic liquid at the titanate interfaces is described 
by an effective one-band Hubbard model on a lattice with N sites which cor- 
responds to the interface Ti02 layer. Each site i identifies a doubled \/2 x y/2 
Ti02 unit cell with two Ti atoms j = 1,2 and four nearest neighbouring oxy- 
gen atoms I — 1,...4 (see Figure [lj, where the cell doubling is introduced 
for the studies of magnetically ordered states on two Ti sublattices. In this 
configuration, the oxygen vacancy corresponds to the elimination of one of 
the oxygens in the Ti 2 04-plaquette. The local disorder induced by an oxygen 
vacancy on the lattice site (i, I) is introduced through the local random fields 
hia,i = her, which shift the electronic 3d xy states of the neighbouring Ti atoms: 



N 




i=l 



tr;j=l,2 




(1) 



1=1 



...4 
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Fig. 1 Scheme of a doubled unit cell on the M02-plane (M=Ti) of a SrTiOa-layer 
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Here d icr ■ are electron creation operators and ni a ,j — d\ a jdi a j are the oc- 
cupation numbers for the self-doped 3d xy electrons with the local energy Ed 
and chemical potential fi; U is the local Hubbard repulsion and t is the ef- 
fective indirect d — d electronic transfer energy. The binary discrete random 
variable Xu = {0, 1} is zero if the oxygen atom is absent in the oxygen po- 
sition of the unit cell i. The two last terms in H describe the magnetic 
splitting = — hi = — h of the local electronic states of Ti due to an oxygen 
vacancy in each of four possible neighbouring (i, I) positions [I51ll7| . For Ti2, 
the nearest neighbouring oxygen atoms belong also to different unit cells with 
the coordinates R;/ = {R,:, R s : — a^, R; + a. y , R; — a. x + a y }. 

In the ferromagnetic state, the local magnetic moment m = rrij = (ntfj) — 
( n ii,j) is defined through the average orbital occupancies of the majority ver- 
sus the minority local spin states which can be expressed via the thermody- 
namic averages of the corresponding double-time one-particle Green functions 

= -m-mdivA^A'vA 1 ')}) US- Hereof (t) anddj,^') 
are the Heisenberg representations of the operators di a .j and d\, a ■, . To cal- 
culate the average electronic orbital occupancies for the Hamiltonian with 
the random configurational variables xu, we employ the coherent potential 
approximation [2T)1I2T] which allows to express the configurationally averaged 

Green functions (G^.j /(T , (i — i')) c through the effective medium Green func- 
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tions R" At-t'): 

N 

f~ijj' _ Tfjf i \ "* T)ja rpi r>gj' io\ 

l = l 

9=1:2 

where the T matrix {T l g } describes the local scattering on the defect potential. 
In the (k, w)-space, the effective medium sublattice Green functions have the 
following form: 

R» a %u) = R» k, W = — — -i jr -J £ , (4) 

2irD a \ uj — x\ ui — x 2 J 

where x\. 2 = (e\ + ef )/2 ± D CT /2, ef. 2 = e rf - /i + U(n s . 1/2 ) + Ef. 2 , and 
£> CT = \J (ef — £r 2 ) 2 + 4i 2 |zk| 2 with Zk = 1 + exp^A^a^) + exp(~-ik y a y ) + 
exp(ik(a a; — a y )) and (n icr;1 / 2 ) = ( n a-.i/2)- The effective self-energies JCj should 
be determined from the equality of the effective-medium propagators and the 
corresponding configurationally averaged Green functions 

< i v(*-*') = (< i v'(*-0)c, (5) 

which is equivalent to the condition {T l g ) c = in the single-site approximation, 
and in our case leads to the following equations for the determination of Ej : 

(Vj°(l-9°Vi a r 1 )c = 0; {J = 1,2), (6) 

where V{ a = h a (A - x a - Xi_a x ,2 - x i+ayA - Xi-a„+ay,3) - Ef and V 2 la = 
h a (A— xa) — E 2 are the random deviations from the effective medium self 
energies due to the local disorder. The local Green functions 

(sj)- 1 = (G^y 1 - EJ (7) 

are determined from the bare Green functions for the stoichiometric lattice: 

(G^-r 1 = cj - e d + fx -U(n & .j) (8) 

which are calculated after the mean-field decoupling of the local Hubbard 
term in H. We note that the mean-field approach allows to capture the main 
features of the ordered states in the correlated systems and is widely used to 
study orbital physics at correlated interfaces [51122]. 

In the ferromagnetically ordered state, we have = E 2 , Gq X = Gq 2 , and 
the problem is reduced to the solution of the equations ^ and the selfconsis- 
tent equations for the magnetization and chemical potential 

m = (n a ,j) - (n 9t j) 

n = (n aJ ) + (n 9ij ), (9) 

which should be considered for a given electron concentration n. The thermo- 
dynamically averaged occupancies (n„j) in this case are expressed through 
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Fig. 2 Local magnetic moment m versus O-vacancy concentration c for different U and elec- 
tron concentrations n. Here h/t = 0.25, ksT/t = 0.01. In plot(b), the electron concentration 
is fixed to n = 0.4. 



the effective-medium Green functions i?J? CT (k, u) and correspond to the con- 
figurationally averaged electronic occupation numbers. 

In the limit of small concentrations of the oxygen vacancies c = (1 — xu) c , 
we can obtain the following expansion for SJ 

Z](u)=4h (7 -^S j (u) + 0(c 2 ), (10) 
J 1 — c 

where Sj(oj) = 1 + h a /(oj — Ed + M — U{n s -j) — h a ). From (TTU1) we see that 
EJ — > for vanishing vacancy concentration c — > 0. 



3 Results and discussion 

The numerical solutions of the equations @ for different vacancy concentra- 
tions c have been analyzed in the low-temperature range 0.0001 < ksT/t < 
0.01 for the electronic concentration n in the range between 0.1 and 0.4 which 
corresponds to the concentrations 6 • 10 13 cm _2 -2 • 10 14 cm -2 of the interface- 
doped polar charge in the Ti02 layers measured in Hall effect experiments j9j 
|5] and estimated from ab-initio calculations [rMTBlfTT] . 

In the regime of weak and intermediate Hubbard correlation energies U/t < 
5, we obtained a monotonic increase of the local magnetization for higher 
vacancy concentrations c, which can be observed in Fig.^a). The value U jt w 
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2 — 5 corresponds to the correlation energies of Ti 3d electrons used in ab-initio 
DFT+U calculations in [T51ITC1IT7] . Fig. [UJa) shows that a robust magnetic 
state with magnetic moments larger than 0.05 fis can be achieved only for the 
high vacancy concentration c > 0.1, which corresponds to vacancy densities 
above a critical value 6 • 10 13 cm~ 2 . This range of vacancy densities was indeed 
analyzed in the ab-initio calculations which explains a good agreement of the 
results obtained in [TMTMTT] with the experimental measurements. 

Furthermore, the analysis of the magnetic moments calculated in the range 
of strong Hubbard correlations U/t > 6 shows the existence of a transition from 
the weak-magnetization regime to the regime of strong magnetism, indicated 
in Fig. [U[b). In the strong-correlation regime, large values of the local mag- 
netic moments of the order 0.3 /is are stabilized already for small c, which 
can be explained by the intrinsic property of the Hubbard model to stabilize 
the ferromagnetic state in the concentration range of approximately quarter- 
filling: the transition from the paramagnetic to the ferromagnetic state upon 
an increase of U can be identified from the magnetic phase diagram for the two- 
dimensional Hubbard model [22 , where the concentration range 0.2 < c < 0.5 
corresponds to the electronic doping levels considered in this work. Considering 
the large values of the magnetic moments, the situation at the cuprate inter- 
faces corresponds rather to the strong-correlation regime. One can expect that 
the strong-correlation regime can be realized, for example, at oxygen-reduced 
cuprate interfaces, which would provide a possibility for a strong enhancement 
of the ferromagnetic state by oxygen vacancies, a scenario which may lead 
to the formation of ferromagmetic regions coexisting with a superconducting 
background. 

Summary 

We considered magnetic states of correlated oxide interfaces, where effective 
charge self-doping and magnetically ordered states emerge due to the electronic 
and ionic reconstructions. Employing the coherent potential approximation to 
the effective one-band Hubbard model with disorder, we analyzed the effect of 
random oxygen vacancies on the formation of two-dimensional magnetism. We 
find that the random vacancies enhance the ferromagnetically ordered state 
and stabilize a robust magnetization above a critical vacancy concentration of 
about c = 0.1. In the strong-correlated regime, we also obtain a nonmomno- 
tonic increase of the magnetization upon an increase of vacancy concentration 
and observe a substantial increase of the magnetic moments. This enhance- 
ment appears due to the intrinsic property of strong electron correlations to 
stabilize a ferromagnetic state at electron doping levels 0.3 — 0.5 electrons per 
unit cell, typical for the polar-doped interfaces [HUH] - Although a mean-field 
evaluation might well overestimate the tendency towards a ferromagnetic state, 
the inhomogeneous (disordered) states considered here rather support ferro- 
magnetic correlations: disorder reduces the kinetic energy and the impurities, 
being effectively magnetic, present a seed for at least short-range magnetism. 
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